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ABSTRACT 


The prediction of ducted fan engine noise using a boundary integral equation method (BIEM) 
is considered. Governing equations for the BIEM are based on linearized acoustics and describe 
the scattering of incident sound by a thin, finite-length cylindrical duct in the presence of a 
uniform axial inflow. A classical boundary value problem (BVP) is derived that includes an 
axisymmetric, locally reacting liner on the duct interior. Using potential theory, the BVP is recast 
as a system of hypersingular boundary integral equations with subsidiary conditions. We describe 
the integral equation derivation and solution procedure in detail. The development of the 
computationally efficient ducted fan noise prediction program TBIEM3D, which implements the 
BIEM, and its utility in conducting parametric noise reduction studies are discussed. Un lik e 
prediction methods based on spinning mode eigenfunction expansions, the BIEM does not require 
the decomposition of the interior acoustic field into its radial and axial components which, for the 
liner case, avoids the solution of a difficult complex eigenvalue problem. Numerical spectral 
studies are presented to illustrate the nexus between the eigenfunction expansion representation 
and BIEM results. We demonstrate BIEM liner capability by examining radiation patterns for 
several cases of practical interest. 
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1. INTRODUCTION 


Advanced analytical tools for predicting the sound radiated from engine ducts facilitate the 
design of active and passive noise abatement technology. To be useful in design studies, 
prediction tools should be fast, versatile, accurate, valid for a wide range of engineering 
situations, and implementable on mainstream computer systems. The ability to compute any 
portion of the sound field without the need to calculate the entire field is an important attribute in 
this regard. Conventional computational approaches such as F in ite Element Methods, CFD, and 
Computational Aeroacoustics (CAA) lack this property which lim its their usefulness for 
conducting parametric noise studies. On the other hand, boundary integral or boundary element 
prediction methods calculate the acoustic field pointwise allowing the designer to compute noise 
only at acoustically sensitive regions of space. 

In this paper, we present a boundary integral equation method (BIEM) for predicting ducted 
fan engine noise. The BIEM is based on the equations of linearized acoustics with uniform inflow 
and predicts the sound scattered by an infinitesimally thin, finite length cylindrical duct that has 
been irradiated by some simple source process. Boundary conditions on the duct interior allow 
for an axisymmetric, axially segmented, locally reactive liner with circumferentially uniform 
impedance. The liner can be positioned anywhere inside the duct (see figure 1). The special case 
of a hard wall interior was considered in reference [1], 

Simple acoustic sources, such as point or line monopoles and dipoles, are used to generate 
incident sound. Source configurations are composed of N symmetrically spaced 
(circumferentially) line or point sources and are situated on a disc perpendicular to the duct axis 
(see figure 1). If the source strengths are constant in time and the source disk spins with uniform 
angular speed U. then classical spinning modes appear inside the duct. Spinning modes can also 
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be generated by nonrotating sources with time harmonic strengths. By suitable choices of 
monopole and dipole strengths, rotating source configurations can be constructed to simulate the 
loading and thickness components of fan noise. 

In the early 1960’s, Tyler and Sofrin first introduced the concept of applying linear infinite 
duct, spinning mode analysis to the prediction of ducted fan engine noise [2]. In this pioneering 
research, simple propagation and radiation models were applied to single radial modes incident on 
the duct inlet. Reflection and inflow effects were ignored and only hard wall boundary conditions 
were considered. Many investigators have since expanded these concepts to include enhanced 
radiation models, inflow effects, and liner capability. The volume of research on this subject is 
considerable and will not be reviewed in any detail here. Comprehensive surveys can be found in 
references [3-5]. 

The finite duct BIEM presented here provides two fundamental improvements over simple 
infinite duct prediction methods. First, the coupling of radiation and propagation for the infinite 
duct methods requires an estimate of the generalized impedances at the inlet and exhaust planes. 
These quantities are then used as input to some radiation model, such as Rayleigh’s formula for 
sound radiation from a flanged cylinder, to calculate the acoustic farfield. For the finite duct 
BIEM, knowledge of the generalized impedances is not required. Reflections at the duct 
openings are determined implicitly. Furthermore, un lik e Rayleigh’s formula, the BIEM accounts 
for edge diffraction and permits noise predictions in the duct shadow region. Second, infinite duct 
methods require the solution of an eigenvalue problem in which the internal acoustic field is 
written as a linear combination of circumferential, radial, and axial eigenfunctions. For the liner 
case, the eigenvalue problem is difficult to solve due to the appearance of Bessel functions with 
complex argument. A body of research, which is reviewed in reference [3], has been devoted to 
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the solution of this very difficult eigenvalue problem. The BIEM does not involve the solution of 
the eigenvalue problem because only a circumferential decomposition of the acoustic field is 
required. If desired, the radial and axial content of the computed BIEM interior field can be 
obtained using Hankel and Fourier transform techniques. 

Application of the BIEM to ducted fan noise prediction is a four step process which we 
summarize below: 

Step 1) In accordance with the preceding assumptions we derive a classical mixed boundary value 
problem (BVP) for the scattered acoustic pressure.. 

Step 2) Using single and double layer Helmholtz potentials, the BVP is converted to a system of 
hypersingular integral equations for the unknown layer densities. 

Step 3) The system of integral equations is solved numerically by the method of collocation in 
which the layer densities are approximated by finite series of orthogonal polynomials. 

Step 4) The scattered sound field is computed pointwise by numerically integrating the Helmholtz 
potential representation. 

Theoretical details of the BIEM form the bulk of this paper. Much analysis is devoted to the 
integral equation development and solution technique. The authors contend that innovations in 
these areas have yielded enhanced numerical procedures that simplify calculations and lead to 
rapid noise predictions. A ducted fan noise prediction program, TBIEM3D [6], that implements 
the BIEM has been developed. Computational results are presented to illustrate TBIEM3D noise 
prediction capabilities and to compare with eigenfunction expansion prediction methods. 
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2. BOUNDARY VALUE PROBLEM DERIVATION 


We consider the scattering of sound by an infinitesimally thin, finite length cylindrical duct in 
the presence of a uniform axial inflow with Mach number M. The duct is irradiated by incident 
sound produced by a collection of N simple point or line monopoles and/or dipoles. Acoustic 
propagation and radiation are based on the assumption of linearity. The N sources have equal 
strengths, are located on a disk perpendicular to the duct axis, and are arranged symmetrically 
about the disk (refer to figure 1). The source disk is centered on the duct axis and either spins or 
is stationary. Source strengths are chosen so that the incident acoustic field can be written as a 
superposition of time harmonic, circumferential modes. 

In the analysis that follows, all quantities are nondimensional: length by r D , mass by p 0 r D \ and 

time by . 

c 


2.1. GOVERNING DIFFERENTIAL EQUATIONS 

We adopt the point of view that the duct is translating in the axial (+z) direction with uniform 
speed U and initially consider an Eulerian description of the acoustic field. The total acoustic 
pressure in the sound field is split into known incident and unknown scattered parts 

p'(r,y/,z,t) = p'(r,ll/,z,t) + p'(r,W,Z,t). ( 1 ) 

In regions of space and time that contain no scattering surfaces, p' is governed by the 



Acoustic pressure and radial velocity are related through the radial component of the acoustic 
momentum equation 
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( 3 ) 


du ' r i dp ' =Q 
dt dr 

In a frame of reference moving with the duct, the symmetry of the source process is such that 
all dependent acoustic variables can be expressed as linear superpositions of time harmonic 
circumferential modes. For example, the scattered pressure has the form 


Ps( r >V>Z,t)= £p;'(r,z)<? 


n (v 7 \i(kt-nNy) 


and the radial component of total acoustic velocity is written 


i'(r,\j/,Z,t) = y j U:(r,Z)e 


n / 7 \ H kt-nNy ) 


where the stretched, moving axial coordinate Z is given by 


z-Mt 


Incident and total acoustic pressures are written similarly. Modal amplitudes in the BIEM are 
calculated term by term. For notational convenience, we drop the superscript n on the modal 
coefficients and define the circumferential mode number m by the equation m = nN . 

If the source disk rotates with angular speed Q. and the strengths of the N sources are time 

TYl 1O1 V 

independent in the moving frame, then k = — and equations (4a-b) yield the classical 

c 

spinning mode representations. In this case, the incident field can be made to simulate that 
produced by an iV-bladed fan. This is accomplished by manipulating monopole and dipole 
strengths to approximate the loading and thickness components of fan noise. For non-spinning 


sources whose strengths are time harmonic with excitation frequency oj, we have k = and 

c 


the time harmonic factor in (4a-b) can be removed from the summations. 



Rather than work with the customary convected wave equation, we simplify the governing 
equations by defining new dependent variables Q , Q s , Q n and V r by the relations 


Q(r,Z) = P(r,Z)e iKMZ 

(6a) 

Q s {r,Z) = P s {r,Z)e iKMZ 

(6b) 

Q i (r,Z) = P i (r,Z)e iKMZ 

(6c) 

V r (r,Z) = U,.(r,Z)e iKMZ 

(6d) 


Combining (4-6) with (2) yields the two dimensional Helmholtz equation 

Qs=o 


~1 d 

( r-] 

o2 2 

d m 

-U 

r dr 

v dr J 

dZ 2 r 2 


+ K~ 


(7) 


for the m-th modal coefficient of scattered pressure. Using (3-6), the m-th radial component of 
the momentum equation (3) can be written as 


-i—Z R 

e *' V,(r,Z) = £ 
M 


Z 

r -i-Z 

e M ^( r,Z')dZ 


dr 


If the duct is stationary, then (7) and (8) reduce to 


"l d 

( r d] 

~\2 2 

d m 

-U 

r dr 

v dr J 

dZ 2 r 2 


P = 0 


and 


U r (r,Z) = ~(r,Z). 


( 8 ) 


(9) 


( 10 ) 


We note that (7-10) are valid for points not lying on the stretched duct. 


2.2. BOUNDARY CONDITIONS 

To meet noise certification requirements, it is necessary to treat the engine duct with noise 
suppression devices. We consider an axisymmetric, locally reactive liner on the duct interior wall 
whose surface admittance is circumferentially uniform. The admittance in the axial direction is 
assumed to be piecewise constant. This property implies that any portion of the duct interior may 
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be lined or hard and allows for inclusion of an axially segmented liner. The hardwall boundary 
condition is used for the duct exterior wall. 

We introduce the concept of a surface function to facilitate discussions of the boundary 
condition and subsequent integral equation derivations. Let /(r,Z) be an arbitrary field function 
and assume that the stretched, moving axial coordinates of the duct trailing and leading edges are 
a and b respectively. Define the surface functions f ± (Z) by the formulas 

/ ± (Z) = lim/(r,Z) Z e (a, b) . (11) 

r— 

Note that the positive (negative) superscript refers to the duct exterior (interior) surface. 

On the outer duct wall the hardwall boundary condition implies that 

V;(Z) = 0 Z g (a,b) . (12a) 

In reference [7] it is shown that if a represents the circumferentially uniform and axially 

piecewise constant specific acoustic admittance on the duct interior surface, then in the stretched, 

moving frame of reference, the surface modal coefficients of radial velocity and pressure satisfy 

the boundary equation 


-i— Z jM (/ d 

e M VrW + ^j- 
(5 K dZ 


b -i—Z ^ 

e M Q 


(Z) = 0 Z g (a,b) . 


(12b) 


In the absence of flow, the boundary conditions reduce to 

u;(z) = 0 Z e (a,b) 


and 


-U;(Z) + aP (Z) = 0 Z e (a,b). 

The momentum equation (8) is used to e lim inate the radial velocity from (12a-b) yielding 


(13a) 

(13b) 


' dQ} 

V dr j 


(. Z) = 0 Z e (a,b), 


-e 




\ dr j 


(Z) + 


iM~a d 
P 3 k dZ 2 


M 


Q 


(Z) = 0 Z e (a,b). 


(14a) 

(14b) 
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and 


r i—r dQ 


e M —(l,Z')dZ' = 0 
dr 


( 15 ) 


Equations (14a-b) are obtained by differentiating (12a-b) and (15) results by evaluating (12a) at 
the trailing edge. These equations are valid for 0 < M < 1 but (15) is satisfied trivially for M — > 0 
and provides no additional information. Note that for M = 0 and a = 0 (hardwall interior) we 
have the classical Neumann boundary conditions. Also, if we add a second coannular duct, then 
boundary conditions similar to (14-15) also apply to the second duct. 


2.3. FARFIELD RADIATION CONDITION AND EDGE CONDITIONS 


In order to have a uniquely solvable BVP we must constrain the behavior of the acoustic 
pressure in the farfield and at the duct edges. To ensure continuity of velocity at the trailing edge 
we impose the Kutta condition 

lim [Q + (Z)-Q-(Z)] = 0. (16) 

Z—ta 1 - J 


For physically reasonable solutions to exist we also require the acoustic pressure to be integrable 
everywhere and particularly in any neighborhood about the leading edge. 

The boundary value problem is finalized by applying the Sommerfeld farfield radiation 
condition, yielding the constraint 



Vp 


^ +iK Q 

dp 


= 0 . 


(17) 


2.4. BVP SUMMARY 


We summarize the above results by listing a complete two-dimensional BVP for the unknown 
scattered acoustic pressure in terms of the known incident pressure. It is assumed that the 
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functions Q , and its derivatives are known, satisfy (17), and are continuous across duct surface. 
For M > 0 


"1 d 

( r _d] 

~\2 2 

0 m 

-U 

r dr 

v dr J 

dZ 1 r 2 


+ K 


Q s = 0 (r, Z) not on duct surface 


Q{r,Z) = Q i (r,Z) + Q,(r,Z) V(r,Z) 

( O/O Y 


dQ 


ydr j 


(Z) = 0 Z e (a,b) 


-e 


■ i s z ('£e'' 

< dr j 


(Z) + 


iM a d 


P 3 k dZ 2 


-i—Z 

e M Q 




(Z) = 0 Z G (a,b) 


r 


e " Z ^-(l,Z')dZ' = 0 
or 


lim i + [e;(z)-e;(z)] = o 

J |<2 S | <791 < 00 91 any neighborhood about leading edge 


lim, Vp| 

p=V r 2 +Z 2 ->oo 


^<2. A 


<?P 




0 . 


(18a) 

(18b) 

(18c) 

(18d) 

(18e) 

(18f) 

(18g) 

(18h) 


§ 3. THE BOUNDARY INTEGRAL EQUATION FORMULATION 
In this section, we convert the two-dimensional boundary value problem (18a-h) to a system 
of boundary integral equations using Helmholtz potentials [8]. This approach is motivated by its 
computational simplicity relative to purely numerical methods based on finite differences or finite 
elements. 

These latter methods involve the calculation and storage of the entire acoustic field and must 
resort to special treatment at the farfield boundary. Consequently, high frequency, farfield 
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predictions require excessive computer memory and computational time. It is shown in reference 
[9] that these drawbacks can be diminished somewhat by coupling nearfield finite-element results 
to Kirchhoff s radiation formula. 

Potential methods, as applied here, involve the pointwise calculation of the acoustic field by 
evaluating duct surface integrals. Noise predictions are calculated at user defined locations only 
and computer storage is min im al. The potential representation satisfies the Sommerfeld radiation 
condition implicitly. Thus, numerical treatment at a fictitious farfield boundary is avoided. Input 
to the potential integrals is obtained by solving a system of hypersingular integral equations for 
unknown surface functions. The d if ficulties associated with solving singular integral equations are 
mitigated by employing innovative analytical and numerical techniques. 

There is extensive research in the literature on the application of integral equation techniques 
to acoustic scattering problems, much of which is discussed in reference [10]. Of particular 
relevance to the BIEM is the work in references [11-14]. The theory and integral equation 
terminology presented in reference [11] provide the foundation for the BIEM theoretical 
discussions. Martinez in references [12-14] has applied singular integral equation methods to the 
study of sound scattered by a thin, finite-length duct. We expand upon this original research by 
applying advanced integral equation solution methods and include the ability to treat arbitrary 
circumferential modes, uniform inflow, enhanced liner models, and a wider range of operating 
parameters. 

3.1. HELMHOLTZ POTENTIAL REPRESENTATION 

The analysis begins with an analytical expression for the Green’s function, G(r,r',Z - Z'), for 
the two-dimensional Helmholtz operator (18a). G satisfies the radiation condition (18h) and can 
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be written 


where 


1 f e~ iKR 

G{r,r',Z - Z') = — cos my r dyr 

2n J R 

0 

R = -Jr 2 +r' 2 —2rr' cosy/ + (Z — Z')~ . 


(19) 

( 20 ) 


To facilitate integral equation analysis and numerical computation, we introduce an integral 
operator notation. For an arbitrary integrable surface function/, define single and double layer 
field operators, S and d, by the equations 


and 


S [f]{r,Z) = jf{Z')s(r,Z-Z')dZ' 


d [fir, Z ) = j f(Z')d(r , Z - Z')dZ ' , 


( 21 ) 

( 22 ) 


where the kernels s and d are given by 

s(r,Z-Z') = G(r,l,Z-Z') (23) 

and 

d(r,Z-Z') = -^(rJ,Z-Z'). (24) 

dr 

The integrals in (21-22) are well defined for points (r,Z) not on the stretched duct surface and 

are equivalent to the classical single and double layer potentials. Care must be exercised when 
evaluating (21-22) on or near the duct surface. 

We define additional field operators by calculating the first radial derivative of S and d. 
Denote by S r and d f the operators 

s r [/](a) = |-{s[/](r,Z)} = J f{Z')s r (r,Z-Z')dZ' (25) 

dr 

and 

d r [/](r,Z) = ~{d[/](r ,Z)} = b jf(Z')d r (r,Z-Z')dZ' , (26) 

dr 

where 
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s r (r,Z-Z') 


ds 

dr 


(r,Z — Z') 


(27) 


and 

d r (r,Z-Z') = ^(r,Z-Z'). (28) 

or 

By applying the Helmholtz potential representation, the scattered pressure can be written as a 
sum of single and double layer Helmholtz potentials yielding 

Qs(r,Z) = s[r /; ](r,Z) + d[r /2 ](r,Z) . (29) 

Equation (29) is a solution of (18a) and (18h). The unknown layer densities q 2 and q 2 are related 
to the jumps in acoustic pressure and its radial derivative across the duct surface. The densities 
are a solution to a system of hypersingular integral equations, which we derive below, and once 
determined provide input to (29) for the pointwise calculation of the acoustic pressure field. 


3.2. SURFACE OPERATOR NOTATION 

To apply the boundary conditions (18c-e) to the potential representation (29), it is necessary 
to evaluate directly the single and double layer potentials and their derivatives on the stretched 
duct surface. The resulting integrals yield one-dimensional surface operators with singular kernels 
that require analytical treatment for their evaluation. 

For Z € \ci,b\ and sufficiently smooth/, we define the surface operators S and D by directly 
evaluating the field operators (21-22) on the duct surface. Thus 


b 


S[/](Z) = j f(Z')S(Z - Z')dZ' 

a 

(30) 

D[/](Z) = J f(Z')D(Z - Z')clZ ' , 

(31) 


a 


where the kernels S and D are given by 

S(Z-Z') = s(l,Z-Z') (32) 
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and 


D(Z-Z') = d(l,Z-Z'). (33) 

Similarly, we form the operators S r and D r by direct evaluation of the differentiated field 
operators (25-26) which yields 


and 


where 


and 


S r [f](Z) = jf(Z')S r (Z-Z')dZ' 

a 

D r [fIZ) = ]f(Z')D,.(Z-Z')dZ' , 

a 

5,(Z-Z') = |i( 1.Z-Z') 

or 

D r (Z-Z') = ^-(l,Z-Z'). 


(34) 

(35) 

(36) 

(37) 


3.3. SURFACE KERNEL PROPERTIES 


The kernels (32-33) and (36-37) are singular for Z-Z' = 0 . It is important both theoretically 
and computationally to determine the precise nature of the singularities. We analytically separate 
the singular portions of the kernels from the bounded parts by applying small argument analysis. 
Details are given in appendix A. In equations (38-41) below, the kernels are written as a sum of 
singular and bounded terms. The bounded parts of the kernels are denoted by a superscript B and 
are written out explicitly in appendix A. 


S(Z-Z') = -^-ln\Z-Z'\ + S B (Z- Z ') 


D(Z -Z') = --j- ln\Z -Z'\ + D b (Z- Z ') 


S,.(Z-Z') = -D(Z-Z') 


(38) 

(39) 

(40) 
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D,(Z-Z') 


1 


(41) 


4(k 2 —m 2 ) + 3 

- + — ln\Z - Z'\ + D? (Z - Z ') 

2k(Z-Z'Y 16ti 

The leading behavior of the kernels (38-40) is logarithmic. Therefore, the associated 
operators are weakly singular. The leading term of the kernel (41) is of the strongly singular 
Hadamard type [15]. Consequently, the integral in (35) is divergent and must be interpreted in 
the Hadamard finite part sense. 

3.4. LAYER CONTINUITY PROPERTIES 

Using the above surface operator notation, we state well-known continuity properties for the 
single and double layer field operators (21-22) and their radial derivatives (25-26) as field points 
approach the stretched duct surface from the interior or exterior of the duct. For a sufficiently 
smooth surface function/, we have the following results [8]: 


lims[/](c,Z) = S[/](Z) 

r— >1 

(42) 

limd[/](r.Z) = +i/(Z) + D[/](Z) 

r-U 2 

(43) 

limS r [/](,-.Z) = +\f(Z) - D[/](Z) 

(44) 

limd r [/](r.Z) = D r [/](Z) 

(45) 


Applying (42-45) to (29) establishes the formulas 

e;(z)-e;(z) = -« 3 (z) (46) 

and 

v dr J 

which relate the layer densities to the jumps in scattered pressure and its radial derivative across 
the stretched duct surface. 


(Z) = -,|(Z) 


(47) 


'dOZ 

dr 


( 2 ) 


17 



3.5. BOUNDARY INTEGRAL EQUATIONS 

By applying the continuity results (42-47), we replace the boundary conditions (18c-e) with an 
equivalent set of integral equations for the unknown surface functions and q , in terms of the 
known incident pressure and its radial derivative. 

The exterior boundary condition (18c) can be written 

-jI + D^,Xz) + D r [ ?2 XZ) = --^-(7,Z) Z e (a,b) , (48) 

where I is the identity operator. This integral equation is hypersingular due to the presence of the 
operator D r . It is worth mentioning here that in many applications, such as lifting wing 
problems, it is customary to work with (48) after performing an axial integration (recall the 
comments following equations (14-15)). In this case, the dominant kernel has a less singular 
Cauchy term. However, no advantage is obtained by this extra calculation because the singular 
operators are computed analytically. Furthermore, for numerical work, we find that working with 
the integrated version of (48) is much more cumbersome. 

The interior boundary condition (18d) is replaced by the sum of (18c) and (18d) yielding the 
system of equations 

(Z) + ^-TX"« Z e -'l(z) = 0 Ze(a,f>) (49) 

fj K (xZj l 
and 

e-(Z) = aU.Z) + %,](Z) + fA + D)[ 9 2 ](Z) Ze(a,b). (50) 

Equations (49-50) can be combined to yield a single equation by computing two tangential 
derivatives of Q . The resulting integral equation contains both Cauchy and Hadamard type 
kernels. It is computationally simpler to introduce a third unknown function, namely Q, and 
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work with the integro-differential equations (49-50). We expound on this in the integral equation 
solution section 4.4. For small Mach number and large k (49-50) can be approximated by the 
much simpler equation 

(/ 3 3 I + iaKS)[q 1 ](Z) + iaK -I + D \[q 7 ](Z) = -iCCKQ^LZ) Z e (a,b) . (51) 

v2 J 

The boundary integral equations are supplemented by (18e), which after substitution of the 
radial derivative of (29) produces the auxiliary equation 

e ^ Z Z,) + S r h> l 1 ' Z ') + d r hi p, Z'^dZ' = 0. (52) 

As mentioned earlier, this equation is satisfied trivially for M = 0 and provides no additional 
information. 

If the interior wall is hard (a = 0), then equations (48-52) reduce to the familiar Neumann 
problem 

q,(Z) = 0 Z e (a,b) , (53) 

D r [q 2 ](Z) = -^(lZ) Z e (a,b) , (54) 

with auxiliary condition 

e ~^ Z l^(lZ') + d r [q 2 ll,Z')\dZ' =0. (55) 

The solution of (53-55) was considered in reference [1]. 

3.6. EDGE CONDITIONS 

The system of integral equations (48-52) admits infinitely many solutions and must be 
supplemented with additional conditions for uniqueness. The as yet unsatisfied BVP edge 
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conditions (18f-g) provide the required constraints by restricting the behavior of the double layer 
density, q 2 , at the duct leading and trailing edges. We omit the edge condition derivations, which 
involves asymptotic expansions of the hypersingular portion of the singular integral equation (48). 
See reference [16] for a detailed discussion of this analysis. 

We finalize the boundary integral equation formulation by stating the asymptotic behavior of 
q 2 near the duct edges: 

q 2 (Z) = 0(VZ - a'j Z —> a + (56a) 

and 

q 2 (Z) = 0\ —^— ) + 0Ub-Z) Z^b~. (56b) 

wfc ~z J 

Note that for M > 0, q 2 has an infinite singularity at the duct leading edge, which is a well-known 
result in lifting surface theory. This term disappears for the no inflow case. 

3.7. BOUNDARY INTEGRAL EQUATION SUMMARY 

We have succeeded in replacing the BVP (18a-h) by an equivalent boundary integral equation 
formulation which we now summarize. 

For 0 < M < 1, the m-th circumferential mode of total acoustic pressure is given by 

2, (r,Z) + s[</ 7 ](r,Z) + d[g 2 ](r,Z) (r,Z) not on stretched duct surface 
Q(r,Z) = < , (57) 

Q i (r,Z) + S[q 1 ](Z) + (+^-l + 2 ](Z) Ze(a,b ) r -> ±1 

V 2 J 

where q 2 and q 2 satisfy 

1) the system of strongly singular integral equations 

-^I + Dk](Z) + D r [ ? ,](Z) = -^(l,Z) Z g {cub) , (58a) 

V 2 J dr 
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-e~' uZ qi (Z) + f e~' uZ Q-\z) = 0 Ze(a,b ) (58b) 

fj K dZ l 

e-(Z) = <2,(i,Z) + % 7 ](Z) + ^I + D^ 2 ](Z) Ze(a,b ) , (58c) 

2) the auxiliary condition 

e ^ Z {^ lZ ^ + S rbJ lZ ') + d rh2](lZ')^dZ' = 0 , (58d) 

and 3) the edge conditions 

q 2 {Z ) = 0{yjZ — a) Z — > a + (58e) 

q 2 (Z) = o[ ^ Ad — \ + QUb-Z) Z-*b~ . (58f) 

\\b Z J 

For coannular ducts, additional integral operators and layer densities must be defined to 
accommodate both ducts. In this case, (57) and (58a-f) have to be coupled to similar sets of 
equations to account for scattering by the second duct. Theoretically, this procedure does not 
change the singular character of the integral equation system nor the solution method. However, 
the computational complexity is nearly double. 

§ 4. INTEGRAL EQUATION SOLUTION 

In this section we present details of the numerical solution of (58a-f). Existence and 
uniqueness results are briefly discussed for the hard wall interior case. A theoretical discussion on 
the solvability of the entire system (58a-f) is lengthy and will be the subject of a future paper by 
the authors. 

Due to edge conditions for the double layer density (58e-f) we are motivated to seek a 
solution to (58a-f) in the form 
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(59) 


q 2 {Z) = A 


Z - a 
\b-Z + 


yj(b — Z)(Z — a) q(Z) 


Z e ( a,b ) . 


In (59), A is an unknown constant and q an unknown Holder continuous function. 

There are no restrictions on the edge behavior of the single layer density. We note however, 
that since the interior wall admittance can be discontinuous it follows from (58b) that q } is also 
discontinuous at the points of discontinuity in a . 

Our solution methodology is best understood by first examining the simple case of a hardwall 
interior with no inflow. We next discuss the hardwall interior case in the presence of inflow, 
which complicates the solution process due to the infinite singularity in the double layer density at 
the duct leading edge. Finally, we expand the method to include the liner case. 


4.1. HARDWALL INTERIOR WITH NO INFLOW 


We simplify the analysis by the introduction of additional operator notation. Let / be an 
arbitrary Holder continuous function and let K denote any of the previously defined integral 
operators, then we define new operators K - by the formulas 


K±[/](Z) = K {b-Zf^(Z'-a)^f (Z) Ze(a,b) 


(60) 


For M = 0 and a = 0, application of (59-60) to (58a-f) produces the single integral equation 

D + r \q](Z) = -^(lZ) Z e (a,b) . (61) 

Equation (61) is a one-dimensional, first kind integral equation with a Hadamard kernel. Integral 
equations of this type appear frequently in the literature and have been applied to problems in 
acoustics, aerodynamics, elasticity, and electrodynamics (see references [17-18], for example). 

Due to the appearance of the weight function *J(b - Z)(Z - a ) , it is natural to seek a solution 
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to (61) in the form 

^4 TT (2Z-a — b\ 

q{ z ) = Y* a ” u » — r • ( 6 2 ) 

It is shown in reference [19] that if the source term in (61) is expandable in a series of second kind 
Chebyshev polynomials, then (61) is uniquely satisfied by a series solution such as (62) except 
possibly at a countable set of eigenfrequencies. In section 5.2, we demonstrate numerically that 
some of the eigenfrequencies of D~ correspond approximately to the transverse resonant 
frequencies of the infinite duct problem. 

In practice, (62) is truncated and the unknown constants a n must be determined numerically, 
usually by some projection technique such as the method of collocation or Galerkin. The 
collocation method is simple to implement and requires fewer calculations than other popular 
methods. For these reasons it is the solution method of choice here. 

To apply the collocation method, we consider the approximate solution 

j(Z)T..l/.f 2Z ~' l " t ) (63) 

V b-a J 

n-0 

and define N? + 1 distinct collocation points 

{Zj} N £ Z^ia.b) Vj . (64) 

We next require the approximate solution to satisfy (61) at the collocation points, which yields the 
linear system of equations 

X/.D IMZ,)=-|(U,) j = l,...,N 2 +l (65) 

n—0 

for the a„’s. Golberg has shown that if the collocation points are chosen as the zeroes of U N +1 , 
then (65) is solvable, except at eigenfrequencies, and the series (63) converges to the actual 
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solution <7 [19]. 

The evaluation of (65) involves the calculation of integrals with finite-part and logarithmic 
integrands. To avoid time consuming numerical integration of these terms, we apply the 
analytical results of appendixes A and B. The bounded portion of the integral operator is 
computed using Gaussian quadrature with weights and nodes based on the second kind 
Chebyshev polynomials. The above combination of collocation points and numerical integration 
scheme yields results as accurate as those obtained by the more computationally intensive 
Galerkin method [19]. 

4.2. HARDWALL INTERIOR WITH INFLOW 

For a = 0 and 0 < M < 1 equations (58a-f) reduce to 

AD r -[i](Z) + D r + M(Z) = -^(i,Z) Z e (a,b) (66) 

or 

and 

a 

e"^ Z (1, z') + Ad; [1](1, Z') + d + r [q](l, Z')]dZ' = 0 . (67) 

In (66-67), A and q(Z) are unknown. We briefly show, under certain conditions, that if the no 
inflow problem (61) is uniquely solvable, then so is (66-67). Furthermore, under the same 
conditions A and q(Z) can be determined separately which we also prove below. 

Define the linear functionals and the known constant C 0 by the relations 

a . K , 

gd [/]= je~ " Z d± [/](;, Z>/Z' (68) 

and 

f -'- z ' dQ 

C 0 =- e M l,Z')dZ ' . (69) 

J dr 

Using (68-69) we rewrite (67) as 
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AiAl] + i*Aq\=C 0 . (70) 

Assume that conditions are such that (61) is uniquely solvable, then the operator has a left 


inverse, which we denote by ^D r J . Consequently, we “solve” (66) for q to yield 

^(Z) = 0j(Z) + A0 2 (Z) Ze(a,b) , 


where 


<t>i(Z) = - D 


-f D rT 7 

\ dQ, 1 

1 r ) 

dr 


(Z) Zs(a,b) , 


and 


02 (Z) = -Dj[i](Z) Ze(a,6) , 

D,=K)X . 


(71) 

(72) 

(73) 

(74) 


Note that and (f) 2 are known functions. 

We obtain an explicit expression for the unknown constant A in terms of known quantities by 
applying (71-74) to (70). This substitution gives 

A = • (75) 

8dl7]-g cM 

If g|([/] ^ [</> 2 ], it follows that A is finite and (71) defines a meaningful expression for q. 

The above results suggest the following algorithm for solving the hardwall interior integral 
equations with inflow: 

Step 1) Solve the no inflow case, i.e., invert D^. 

Step 2) Calculate A via (72-75). 

Step 3) Calculate q using (71) and (75). 

Theoretically, this approach is equivalent to the no inflow case, which implies that the results 
of Golberg apply [19]. Furthermore, since the computational time required for the calculation of 
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the functionals (68) is negligible compared to the time needed to compute the matrix for D*, it 

follows that the above algorithm is also computationally equivalent to the no inflow case. 

In reference [1], a slightly different approach was used to separate A and q. The resulting 
algorithms were subsequently shown to be convergent [20]. 

4.3. LINED INTERIOR WITH NO INFLOW 

In this case, (58a-f) reduce to the system of integral equations 

( 1 \ 

(I + ia kS)[< 7 ; ](Z) + ia k - 1 + + D + \q]{Z) = -ia kQ { (1, Z) Ze (a, b ) . (76a) 

\2 J 

-jI + d)[ ? ,](Z) + D;MZ) = -^(;,Z) Ze(a,b) (76b) 

We defer a detailed theoretical discussion of the singular system (76a-b) to a future publication 
and, instead, focus on the solution technique. 

To solve (76a-b) numerically we again apply the collocation method. The double layer 
density q is approximated as in the no inflow, hardwall case. Since there are no restrictions on the 
single layer edge behavior and the weight function for the Legendre polynomials is unity we seek 
an approximation to q ; of the form 

,-(Z)Tw(^) • (77) 

V b-a J 

n=0 

Further motivation is provided by the fact that the numerical accuracy achievable by the expansion 
(77) is commensurate with that for the double layer density approximation (63). 

We collocate (76a) at the zeroes of P Ni+1 , [T ; ] , and (76b) at the zeroes of U Ni+I , 

|Z ; . } . ; , yielding the IN : N 2 ) x (NjN 2 ) linear system 
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(78a) 


Nj Nj 

K (I + iaicS)[P„ ](Yj ) + iccK ^ «„ 

n=0 n=0 

= -iccKQ i (l,Y j ) j = 



jI ++D+ ][f/„](F 7 ) 


“£ b * [ 2 1 + D l Z J ) ' + £ ' fl » D r ) ; = - ^ Z j) 


j = 2 (78b) 


for the unknown expansion coefficients. No theoretical results exist regarding the solvability and 
convergence of this numerical method. However, computational evidence suggests that the linear 
system (78a-b) is stable and the discretized solutions converge to the actual solution of the system 
of integral equations (76a-b). 


4.4. LINED INTERIOR WITH INFLOW 
Define the functional g s by the formula 

a K 

g s [/]= J e~^ Z S r [f]{l,Z')dZ' . (79) 

For a lined interior and 0 < M < 1 , equations (58a-f) assume their most complex form: 

-e"« Z qi (Z) + f e"« Z Q-\z) = 0 Ze(a,b) (80a) 

p K dZ l 

Q(Z) = Q i (l,Z) + S[q 1 ](Z) + 

f / x f 1 \ (80b) 

A[-I“ + D J [l](Z) + )-I + + D + pjz) Z s (a,b) 

- |l + Dj[ 9l ](Z) + D+[ 9 ](Z)+AD i :[l](Z) = -^(l,Z) Ze(a,b) (80c) 

S s [qi] + A ga[l] + gM= C o (80d) 

Equations (80a-b) can be combined to form one equation by substituting (80a) into (80b). 
The resulting hypersingular integral equation contains two axial derivatives of the single and 
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double layer surface operators. It is far simpler computationally to introduce a third unknown 
function Q , which we again approximate by a finite sum of Legendre polynomials. A linear 
system is obtained by applying the usual collocation method and the unknown constant A is 
determined by a manner similar to that used for the hardwall case. 

We simplify the solution process even further by replacing (58b-c) with the approximate 
equation (51). Thus, (80a-b) reduce to 


(j /8 J I + /a/cS)[qj](Z) + ia/<A — I +D Y/](Z) + 

v 2 J 

ia trf 1 1 + + D + \q](Z) = -ia kQ (i, Z) Ze ( a , b ) 


(81) 


This equation is valid for small values of the parameter and is similar to the no inflow 

integral equation (76a). The system of equations (81) and (80c-d) is then solved by first 
separating the determination of the constant A and then applying the collocation scheme described 
in the previous section. 


§ 5. RESULTS 

In this section, computational results are presented that demonstrate BIEM features. We first 
discuss the ducted fan noise prediction code TBIEM3D [6] and its application to parametric noise 
reduction studies. Next, the relationship between the finite duct BIEM and infinite duct 
prediction methods is established by analyzing the spectral content of the BIEM interior pressure 
field. Lastly, we illustrate BIEM liner capability by conducting a numerical experiment in which 
the effects of interior admittance distributions on radiated sound are studied. 
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5.1. THE DUCTED FAN NOISE PREDICTION CODE TBIEM3D 


A PC based computer program, TBIEM3D (Thin duct, BIEM, Three-Dimensional), has been 
developed which incorporates the BIEM solution features [6]. Results from a typical TBIEM3D 
calculation are shown in figure 2. In this simulation, 20 spinning point axial dipoles, situated on a 
disk in the middle of a short duct, generate the twentieth circumferential mode of acoustic 
pressure and its harmonics. The duct walls are hard and sound propagates in the presence of a 
uniform axial flow with M = 0.4. This configuration was chosen to approximate the thrust 
component of loading noise produced by a 20-bladed fan. Figure 2a shows a cut-a-way view of 
the engine duct in which contours of instantaneous acoustic pressure for the m = 20 
circumferential mode are plotted in the duct interior. Conditions are such that only one radial 
mode is cut-on. Doppler effects and reflections from the duct openings are included. Portions of 
sound reaching the duct ends radiate to the farfield as indicated in figure 2b. The combined 
acoustic field in figure 2 is composed of 10,000 field calculations and consumed approximately 4- 
5 minutes computational time on a laptop personal computer with a Pentium 133 processor and 
32 megabytes of random access memory. 

To demonstrate the effectiveness of TBIEM3D as a noise prediction tool, the authors and 
coworkers have applied the code in several engineering studies. In reference [21], TBIEM3D 
was used to quantify the concept of generalized radiation impedance and to assess the accuracy of 
combining infinite duct propagation results with Rayleigh’s radiation formula. Results from this 
research clearly illustrate the similarities and differences between the BIEM and infinite duct 
prediction methods. 

TBIEM3D has been applied to both passive and active noise control research [22-25]. For 
given liner length and location, optimal liner impedances were calculated for a host of operating 
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situations in which circumferential mode number, inflow Mach number, and excitation frequency 
were varied [22]. Optimal liner impedance is also a function of the designer’s definition of noise 
attenuation. In this regard, TBIEM3D can be easily adapted to accommodate noise metrics that 
are based on the radiated or internal acoustic pressure fields. For example, in reference [22] liner 
impedance was chosen so that the max im um sound pressure level along an arc of farfield 
observers in the direction of the principal radial mode was min im ized. Active noise control 
research based on TBIEM3D has focused on using simple point control sources and a 
feedforward algorithm to min im ize the sound along a particular farfield direction produced by a 
given source configuration [23]. 

A two dimensional version of the code has been used to examine the effect of a scarf inlet on 
farfield noise radiation from a short duct [24-25]. Results from these studies indicate that the 
scarf inlet is very effective at redirecting high order transverse modes away from noise sensitive 
regions over a wide range of excitation frequencies. This phenomenon is well known at high 
frequency [26]. Calculations suggest that the scarf inlet also redirects lower order transverse 
modes but to a lesser degree. In the absence of flow, the scarf inlet has little affect on aft 
radiation patterns. However, the amount of distortion in the aft region due to the scarf inlet 
appears to increase with increasing inflow Mach number. It was further shown that the interior 
acoustic field may be seriously altered by the presence of a scarf inlet, which provides an 
additional mechanism for passive noise control. 

5.2. HARDWALL INTERIOR: SPECTRAL ANALYSIS 

In this two part numerical experiment, we analyze the spectral contents of the BIEM hardwall 
interior field and relate the findings to the usual infinite duct eigenfunction expansion. To 
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simulate infinite duct conditions with TBIEM3D, we consider a long duct (L d = 6.0 ) with the 
source disk located at one end. 

For an infinitely long hardwall duct with unit radius, it is well-known that, due to the 
phenomenon of transverse resonance, the interior Helmholtz problem is unsolvable at the cut-on 

wavenumbers \k r 1 , where J' m (k '” n ) = 0. It is intuitive to assume that the discretized linear 

L J n=] m 

system for the long, finite duct problem should experience some ill-conditioning near the infinite 
duct eigenwavenumbers. 

In figure 3, we plot the L t condition number of the D+ matrix (see equation 65) as a function 
of excitation wavenumber, k, for the circumferential mode numbers m = 4 and m = 10. The 
location of the infinite duct resonant wavenumbers are denoted by dashed lines. As expected, 
serious ill-conditioning occurs near the cut-on wavenumbers. Also present are lesser spikes of ill- 
conditioning. The cause of the finite duct ill-conditioning is probably due to the organ-pipe 
resonance phenomenon. This subject and the viability of the BIEM solution at these organ-pipe 
eigenwavenumbers will be addressed in a future paper by the authors. 

Part two of this experiment involves comparing BIEM interior axial wavenumber spectra with 
the discrete infinite duct results. The infinite duct axial wavenumbers are given by the formula [3] 



The duct mode (m,n) is cut-on if k > k . To compute the axial spectrum of the BIEM interior 
field, the acoustic pressure is calculated along a line parallel to the duct axis and analyzed using 
FFT techniques. Axial BIEM spectra for two separate cases are presented in figure 4 and 
compared to the in fin ite duct spectra. Operating conditions were chosen so that three radial 
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modes were cut-on. There are two important observations from these results. First, the locations 
of the BIEM spectral peaks are in excellent agreement with the infinite duct axial wavenumbers. 
Second, reflections at the duct ends are included and the relative strengths of propagated and 
reflected modes can be easily determined using BIEM results. 

5.3. SAMPLE LINER CALCULATIONS 

We illustrate BIEM liner capabilities by studying the effect of interior admittance distribution 
on noise radiation patterns. For this virtual experiment, we consider a short duct, L D = 0.5, with 
the source disk located in the middle of the duct. Twenty spinning point axial dipoles generate 
incident sound in the presence of a uniform inflow with M = 0.4. The rotational speed of the 
source disk is chosen such that one radial mode is cut-on. 

In figure 5, contours of sound pressure level are plotted in the acoustic field for four different 
admittance distributions. Figure 5a shows the scattered field for a hardwall interior. Three 
different liner configurations are considered in figures 5b-d. As shown in figure 5b, a portion of 
the inlet is lined ( L = 0.2 ) with admittance chosen to be optimal in the sense of reference [22]. 
Significant noise reduction in the forward region is attained. Figure 5c shows the results of lining 
a portion of the exhaust duct. For this case, the admittance is not optimal. Some noise reduction 
in the aft region is apparent with little or no reduction in the forward region. In figure 5d, both 
the inlet and exhaust are lined resulting in substantial noise reduction both forward and aft. Each 
case presented here required several minutes computational time on a Pentium 133 personal 
computer. 

§ 6. CONCLUDING REMARKS 

A boundary integral equation method for the prediction of ducted fan engine noise has been 
presented. Based on the equations of linearized acoustics, a classical boundary value problem was 
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derived in its entirety. The ability to treat uniform, axial inflow and locally reacting, axially 
segmented liners is included. Using Helmholtz potentials, the BVP was recast in terms of a 
system of hypersingular boundary integral equations. The extension of the BIEM to coannular 
ducts was described. 

Innovations in the theoretical approach to solving the integral equation system and in the 
analysis of the integral equation kernels have produced a computationally efficient solution 
procedure. The ducted fan noise prediction program TBIEM3D incorporates the BIEM solution 
features and rapidly predicts the sound field produced by the scattering of incident sound by a 
finite length cylindrical duct. 

Unlike infinite duct prediction methods, the BIEM does not require the decomposition of the 
acoustic field into its radial and axial components. However, spectral evidence presented in the 
results section clearly establishes the connection between the finite duct BIEM and infinite duct 
prediction methods. That is, for the hardwall case, the BIEM interior field displays the axial and 
radial spectral characteristics that are predicted by the infinite duct theory. Furthermore, 
reflections due to the duct openings are accounted for implicitly without resorting to approximate 
measures. It was also demonstrated that noise predictions involving lined interiors are easily 
calculated without the need to solve the complex eigenvalue problem. 

It has been shown that TBIEM3D is accurate, executable on inexpensive computer platforms, 
and applicable over a wide range of operating conditions. Code versatility, ease-of-use, and rapid 
prediction capability are necessary features for conducting cost-efficient parametric noise 
reduction analyses. In this regard, the effectiveness of TBIEM3D as an engineering analysis tool 
has been amply demonstrated both here and in the literature [21-25]. 
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§ 8. APPENDICES 


8.1. APPENDIX A: SINGULAR KERNEL ANALYSIS 


We derive the singular surface kernel results (38-41) by applying small argument analysis to 
the Green’s function (19) and its derivatives. Since no simple closed form expression exists for 
(19), it must be evaluated numerically. Due to the singularities, conventional numerical 
integration fails when evaluating the Green’s function on the stretched duct surface. Thus, 
analytical methods are required for resolution of the singular terms. For each kernel, the singular 
portion is separated from the bounded part by subtracting the singular terms from the Green’s 
function integrand. The resulting singular integrals are evaluated in closed form and the remaining 
integrands are continuous everywhere and present no numerical problems. 

From (32) and (19), we have 


where 


S(Z-Z') = — 
’ 2k 


-iKR,, 


cos m i i> cl y/ 

R, 


'■a 


R n = 


2 sin — ] +(Z-Z') 2 . 
2 ) 


(Al) 


(A2) 


The integrand in (Al) is O 



{Rn) 


for small y/ and Z-Z’ 


Therefore, we seek a function g s (y/) 


such that 


h s (yr, Z - Z') = cos my/ e~ iKR " -g s (y/) = 0(R n ) yt\Z-Z'\^0 (A3) 


and the integral 


8s(Y) 


R 


dyr 


11 


is known in terms of elementary functions. One such choice is 


(A4) 
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8s(.V) = c os y > 


which yields 


r ¥ 
cos — 

2 

R„ 


Consequently, 


where 


dy/ = — In Z-Z' + In 


1 


2 + J4 + (Z-Z') 2 


S{Z - Z') = ln|Z - Zl + S B (Z - Z') 

2k 


S B (Z - Z') = — 
v ’ 2k 


hJy/,Z-Z') 1 
sKY —d\j/ + — In (2 + IV) 




n 


2k 


and 


W = J4 + (Z-Z'Y 


We apply a similar analysis to the double layer surface kernel to obtain 


D(Z - Z') = — — ln|Z - Zl + D b (Z - Z') 
4k 


where 


and 


D b (Z-Z') = — 
4k 


h D (ys,z- z') 


h n (y/,Z - Z') 1 x 1 

DKY — -dy/ + — In (2+W) 
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4k 


2kW 


f 2 sin — 1 

2 r 
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cosml/r(i + iKRjj)e ,KR " -cos^- 


Analysis for the hypersingular kernel (41) proceeds like before producing 

’H. 


r>,(7-7') = 


1 .4 


2 ) 

K — 111 


2k(Z- Z'Y 16k 

The expressions for D B are lengthy and will be given in two parts: 

D b (Z - Z ') = D b {Z - Z ') + D b (Z - Z') 


ln|Z ~ Z'\ + D b (Z - Z') . 


where 


D‘(Z-Z’) = ~Y 

2k 


h ^ Z - Z X, W - 2 ^ H 2 + W)Z + 2C{2+W) 
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4k 


2kW(2+W) ’ 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 
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2 v ' 8 tt /?, 5 , 871 y ’ 4 kW 3 


(A16) 


hj(y/,Z-Z') = cosm i/a cos y/(l + iKR n )e ,KR " -cos— i + C 2sin— 1 +— i? 2 ; , (A17) 

2 \ 2 J 2 


XI/ i 

h 2 (lj/,Z-Z') = 2 sin— cosmy/(^3 + 3iKR,j - K 2 R 2 j)e~ ,KRn - 3 cos- 


(A18) 


2 + 4hT 


(A19) 


8.2. APPENDIX B: SINGULAR INTEGRAL EVALUATIONS 

In solving the system of integral equations (58a-g), hypersingular and logarithmic integrals 
involving Legendre and second kind Chebyshev polynomials are encountered. We list without 
proof several pertinent integration results, which can be obtained by applying the Plemelj- 
Sokhotski theorem and its logarithmic analog [27]. Some of these results are well-known and 
appear in the literature (see reference [28], for example). 
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n > 0 (B3) 
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§ 9. NOTATION 

9.1. PARAMETERS AND VARIABLES 

Denotes that a quantity is dimensional when appearing over a variable 
a Axial coordinate of duct trailing edge in stretched, moving frame 

b Axial coordinate of duct leading edge in stretched, moving frame 

c Ambient sound speed 

J m First kind Bessel function of order m 

k = -ZZ Characteristic wave number 

c 
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Radial wavenumbers 


7 mn 

K 

k'™ Axial wavenumbers 

L Liner length 

L d Ratio of duct length to duct diameter 

m Circumferential mode number 

M = — Inflow Mach number 

c 

n Radial mode number 

N Number of sources 

p' Eulerian description of total acoustic pressure field 

p' Eulerian description of incident acoustic pressure field 

p' Eulerian description of scattered acoustic pressure field 

P n Legendre polynomials 

Q n Associated Legendre functions of the second kind of order zero 

(r,\j/,z) Cylindrical coordinates in Eulerian frame 
r D Duct radius 

t, t Time 

T n Chebyshev polynomials of the first kind 

u' r Eulerian description of the radial component of acoustic velocity 
U Axial speed of duct 

U n Chebyshev polynomials of the second kind 

Z Axial coordinate in stretched, moving frame 
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a Specific acoustic admittance on duct interior surface 

P = y/l-M 2 Stretching parameter 

K = k/ P Stretched characteristic wave number 

p 0 Ambient density 

cb Excitation frequency (radians/second) 

Pi Shaft speed (radians/second) 

9.2. INTEGRAL EQUATION OPERATORS AND KERNELS 

d, d Double layer field operator and kernel 

D, D Double layer surface operator and kernel 

D b Bounded portion of double layer surface kernel 

d r , d r Radial derivative of double layer field operator and kernel 

D r D r Radial derivative of double layer surface operator and kernel 

l) B Bounded portion of double layer radial derivative surface kernel 

S, 5 Single layer field operator and kernel 

S, S Single layer surface operator and kernel 

S B Bounded portion of single layer surface kernel 

S r , s r Radial derivative of single layer field operator and kernel 

S r , S r Radial derivative of single layer surface operator and kernel 
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Figure 1 : BIEM Geometry - Source Plane Located at z-0 
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Figure 2: Sample TBIEM3D Results (Hardwall Interior) 
Contours of Instantaneous Acoustic Pressure (Pa) 
a) Cut-Away View b) Side View 
N = 20 m = 20 k = 22 M = 0.4 L D = 0.5 
One Radial Mode Cut-on 
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Dashed Lines Indicate Infinite Duct Radial Eigenvalues 


Figure 3: L } Condition Number of D* Matrix as a Function of Wavenumber 

a) m = 4 b) m = 10 
M = 0.0 L d = 6.0 
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Figure 4: Axial Wavenumber Spectra for Hardwall Interior 
a) m = 4 k = 12.0 b)m=10 k = 20.0 

M = 0.4 L d = 6.0 

BIEM Infinite Duct 
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Figure 5: Sound Pressure Levels for Various Liner Admittance Distributions 
a) Hardwall Interior b) Lined Inlet c) Lined Exhaust d) Lined Inlet and Exhaust 

m = 20 k = 24.0 M = 0.4 L D = 0.5 





